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ABSTRACT 

Results from exoplanet surveys indicate that small planets (super-Earth size and below) are 
abundant in our Galaxy. However, little is known about their interiors and atmospheres. There 
is therefore a need to find small planets transiting bright stars, which would enable a detailed 
characterisation of this population of objects. We present the results of a search for the transit 
of the Earth-mass exoplanet Alpha Centauri B b with the Hubble Space Telescope (HST). We 
observed Alpha Centauri B twice in 2013 and 2014 for a total of 40 hours. We achieve a 
precision of 115 ppm per 6-s exposure time in a highly-saturated regime, which is found to 
be consistent across HST orbits. We rule out the transiting nature of Alpha Centauri B b with 
the orbital parameters published in the literature at 96.6% confidence. We find in our data a 
single transit-like event that could be associated to another Earth-size planet in the system, 
on a longer period orbit. Our program demonstrates the ability of HST to obtain consistent, 
high-precision photometry of saturated stars over 26 hours of continuous observations. 
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1 INTRODUCTION 

Ground- and space-based observations have been providing a 
wealth of data on exoplanets covering a wide parameter regime. On 
one hand, dozens of detected hot Jupiters with high-SNR observa¬ 
tions have expanded our knowledge of giant, irradiated exoplan¬ 
ets (e.g. Seager & Deming 2010; Heng & Showman 2014). On the 
other hand lie the smaller, lower mass objects, of which terrestrial 
planets are a subset and about which very little is known. Both Ke¬ 
pler (Batalha et al. 2013) and radial-velocity surveys (Howard et al. 
2010; Mayor et al. 2011) have been detecting an emerging popula¬ 
tion of these small, possibly terrestrial planets that appear to be 


ubiquitous. Very little is known about the interiors of small planets 
due to compositional degeneracies and their atmospheric proper¬ 
ties are challenging to probe. The search for and characterisation 
of terrestrial planets today represent the cutting edge of the exo¬ 
planet field. The harvest of results obtained on transiting giant plan¬ 
ets shows the pathway to improve our knowledge of small planets: 
we have to find terrestrial exoplanets that transit bright stars, for 
which both masses and radii can be precisely measured. Despite 
the large number of Earth-size candidates found by Kepler , most 
of their host stars are too faint to precisely constrain the plane¬ 
tary masses using radial-velocity techniques. For most of the plan¬ 
ets, we know only their sizes, not their masses or densities. How¬ 
ever, radial-velocity surveys target close, bright stars. Thus, pre¬ 
cisely monitoring the expected transit window of radial-velocity 
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discovered planets (Gillon et al. 2011) could potentially result in 
the detection of a transit, the derivation of the planetary radius 
and advanced follow-up techniques. Successful cases are, among 
others, GJ436 b (Gillon et al. 2007), 55 Cnc e (Demory et al. 2011; 
Winn et al. 2011) and HD 97658 b (Dragomir et al. 2013). 

The recent discovery of a Cen B b with the HARPS spectro¬ 
graph (Dumusque et al. 2012) is a significant achievement towards 
the discovery of an Earth twin, by demonstrating the ability of the 
radial-velocity technique to detect Earth-mass companions. With a 
Doppler semi-amplitude of only 0.51 ±0.04 m/s, a Cen Bb is the 
lowest RV-amplitude planet found as of today. Four years of data 
and a careful analysis were necessary for the authors to characterise 
the stellar oscillation modes, rotational-induced and long-term ac¬ 
tivity as well as the binary orbital motion. These critical steps re¬ 
vealed the planetary signal and underlined the need of simultane¬ 
ously fitting the stellar and planetary signals to push the precision 
of radial-velocity surveys (e.g. Boisse et al. 2010). The complexity 
of the analysis, coupled to the weak planetary signal could be seen 
as the main challenge of this major discovery (Hatzes 2013), call¬ 
ing for an independent confirmation. The motivation of the present 
study was to carry out a photometric confirmation of this planet 
with the Hubble Space Telescope (HST). With a geometric prob¬ 
ability of 9.5%, the detection of a transit of a Cen B b would be 
first and foremost a clear confirmation that the planet is real, and 
would represent an opportunity to precisely constrain the density 
of an Earth-mass planet. 


2 OBSERVATIONS AND DATA REDUCTION 

We first observed a Cen B (V = 1.33, B-V = 0.90, K1V) quasi- 
continuously for 16 orbits (26 hours) on July 7-8 2013 with the 
Space Telescope Imaging Spectrograph (STIS). This long moni¬ 
toring sequence allowed us to cover —96% of the transit window 
allowed by the radial-velocity orbital solution of Dumusque et al. 
( 2012 ). 

The brighter a Cen A (V = -0.01, B-V = 0.69, G2 V) was 
4.5” away from the B component in July 2013. The small apparent 
separation and respective brightnesses of both components of the 
a Cen system required tailored observational settings. In particular, 
the stellar brightness forced us to saturate the exposures; not satu¬ 
rating them would have meant setting short exposure times com¬ 
mensurable with the detector shutter closing time, yielding high 
uncertainties on the actual values of the exposure time (shutter¬ 
timing jitter). While previous studies found that saturated spec¬ 
tra do not preclude achieving high signal-to-noise ratio photome¬ 
try on bright stars (Gilliland 1999; Gilliland et al. 1999), saturating 
a Cen B implies that a Cen A would saturate even more, making it 
impossible to have both stars in the spectrograph slit without hav¬ 
ing their blending charges overlapping. To avoid this situation, we 
positioned the spectrograph slit on a Cen B and requested a range 
of telescope roll angles such as a Cen A would be masked out and 
aligned in the dispersion direction. This configuration indeed min¬ 
imises any contamination by a Cen A and the potential impact of 
its 45-degree diffraction spikes. HST, however, does not actually 
roll on request but presents varying roll angle as a function of time. 
In this case, it was possible to schedule observations of a Cen B 
in the HST continuous viewing zone (free of Earth occultations) 
while committing with the selected range of telescope roll angles. 
The target acquisition was first performed on the brighter a Cen A 
before offsetting to a Cen B. The target was then carefully centred 
in the 52" x 0.05" long slit. A wider slit (52" x 2") was finally 



Figure 1. 2D raw extracted spectrum. This spectrum is a sample from our 
program PID13180. The degree of saturation as well as scattered light from 
the star are both clearly visible. 


used to minimise slit losses for science exposures. We show for 
illustration purposes one spectrum of our program in Figure 1. 

We then re-observed a Cen B for 9 orbits on July 27-28 2014 
to confirm a transit-shaped event detected in the July 2013 dataset. 
We provide details about our data reduction and analysis in the fol¬ 
lowing section. 

2.1 July 2013 dataset 

We observed a Cen B using STIS/CCD with the G750M grating 
and an exposure time of 6 s. We use the setting 6094, which has 
a wavelength domain ranging from 5812 to 6380Aand a resolution 
of -5500. 

2.1.1 Data Reduction 

The starting point of our data reduction consists of 11 flat-fielded 
science files (FLT) that are the output of the STIS calibration 
pipeline CALSTIS. These data are available on MAST under pro¬ 
gram ID 13180. Each FLT file is a FITS data cube that contains 
a 2D spectrum for each exposure (2087 in total). All spectra are 
1024 pixels in the dispersion direction and 300 pixels in the cross¬ 
dispersion direction. Since a Cen B saturates the detector in the 
adopted 6-s exposure time, this larger acquisition window allowed 
us to integrate those columns that show bleeding up to —60 pix¬ 
els along the cross-dispersion axis. We identify and discard cos¬ 
mic ray hits and bad pixels on each calibrated frame. The HJD 
mid-exposure time is obtained from the file headers and converted 
to BJDtds using existing routines (Eastman et al. 2010). We use 
custom-built procedures to perform two one-dimensional spectral 
extractions using rectangular apertures with 120 and 180 pixels in 
the cross-dispersion axis centred on the spectrum. We perform the 
extraction with two different apertures initially to ensure that the 
photometry is insensitive to the aperture size. We apply these ex¬ 
tractions to the 11 FLT data cubes. We hold both the aperture size 
and position fixed for all spectral frames. We finally sum the spectra 
to build a high S/N ratio white light-curve for the photometric anal¬ 
ysis. We collect —6 10 9 e~ per 6-s exposure time. We find no ev¬ 
idence of contamination from the A component. Visual inspection 
of the raw photometric time-series shows a —0.4% increase in flux 
for the first 0.6 days combined to a sinusoidal signal that matches 
HST’s orbital period. The raw photometry is shown in Figure 2. 

Tracking the spectral shifts on the detector is rendered diffi¬ 
cult because of the level of saturation. To monitor displacements in 
the dispersion direction, we locate and fit the Na doublet lines 40 
pixels from the spectral intensity centreline in the cross-dispersion 
direction to avoid performing the fitting in the saturated region. In 
the cross-dispersion direction, a section of the spectrum (at a given 
wavelength) can be approximated by a boxcar function, which is 1 
for saturated pixels and 0 elsewhere, that we fit in seven locations 
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Figure 2. Raw HST/STIS photometry. The single-baseline best-fit MCMC model is superimposed in red. The dotted lines indicate the location and width of 
the transit pattern described in Section 2.1.2. 


along the dispersion direction of the spectrum. We use these mul¬ 
tiple determinations to estimate an average value for the spectrum 
centre for each frame. We respectively detect drifts of 0.33 pixels 
RMS in the dispersion direction and 0.11 pixels RMS in the cross¬ 
dispersion direction. We show in Figures 3 and 4 the spectral shifts 
in both directions during the entire visit. 

We measure the median background level in two 10x1022 
pixel regions located at the subarray top edge and in another area 
of the same size at the bottom edge of the spectrum. Time-series 
extracted from both regions have a slightly different shape but their 
Fourier spectra have the same characteristic frequencies. We find 
a background level of about 54 e~ /pixel during the course of our 
observations, which is likely due to scattered light from a Cen B. 
Examination of the background time-series reveals a clear period¬ 
icity at FIST’s orbital period as expected, that accounts for most of 
the measured background RMS at the 1% level. The residual RMS 
for the background time-series is 0.02%. We assess the incidence of 
the background on the time-series by subtracting its median value 
from each pixel in the aperture window. We find no effect on the 
photometry. 

In the next step we focus on the impact of the flat-field correc¬ 
tion on the photometric time-series. Calibrated 2D spectra available 


on MAST 1 have been reduced using the X6417097O-PFL.FITS 
pixel-to-pixel flat file, whose RMS across the full frame is ~ 1 %. As 
previously discussed (Gilliland 1999), the STIS CCD has a very flat 
response across its wavelength range. The flat field noise per pixel 
drift is estimated by dividing its RMS by the square root of the 
total number of pixels over which the incoming flux is spread on 
the detector. The stellar flux is spread over ~60 pixels in the cross¬ 
dispersion direction and 1022 pixels in the dispersion axis (the two 
extreme columns are trimmed). The noise induced by the imperfect 
pixel-to-pixel flat is thus ~38 ppm per pixel shift, which is below 
the measured individual errors of 115 ppm per 6-s exposure. As 
a second check, we also recalibrate the data without the flat-field 
correction, which results in photometric time-series that are indis¬ 
tinguishable from the ones including the pixel-to-pixel sensitivity 
correction. 

We measure the overscan level in a fixed 12x120 pixel re¬ 
gion located at the extreme left of the RAW files (as both overscan 
regions are trimmed during the 2D spectral calibration). Overscan 
time-series exhibit an asymptotic decrease reaching a stable level 
0.6 days after the start of the observations. This trend is combined 


1 http://mast.stsci.edu 
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Figure 3. Spectral shifts in X. Displacement of the spectrum in the dispersion direction as calibrated using the Sodium doublet. The dotted lines indicate the 
location and width of the transit pattern described in Section 2.1.2. 



Figure 4. Spectral shifts in Y. Displacement of the spectrum in the cross-dispersion direction. The dotted lines indicate the location and width of the transit 
pattern described in Section 2.1.2. 


with a higher frequency variation matching HST’s orbital phase 96- 
min duration. 

For each spectrum we record the timestamp, normalised 
flux, photometric error, spectral shift in the dispersion and cross¬ 
dispersion directions, median background and overscan values as 
well as HST’s orbital phase at which each spectrum has been ac¬ 
quired. 


2.1.2 Photometric Analysis 

We use the MCMC algorithm implementation already presented 
in the literature (e.g. Gillon et al. 2012). Inputs to the MCMC are 
the photometric time-series obtained during the data reduction de¬ 
scribed above. Photometric baseline model coefficients used for de¬ 
trending are determined at each step of the MCMC procedure using 
a singular value decomposition method (Press et al. 1992). The re¬ 
sulting coefficients are then used to correct the raw photometric 
lightcurves. 

The baseline model for the full, 26-hr long time-series (see 
Figure 2) consists of a second order logarithmic ramp (e.g. 
Knutson et al. 2008; Demory et al. 2011), combined with a lin¬ 
ear function of time and a fourth-order polynomial of HST s or¬ 
bital phase. The choice of the function best fitting the HST orbital 
phase effect, usually attributed to telescope breathing, is commonly 


based on half HST orbits (e.g. Sing et al. 2009; Huitson et al. 2012; 
Evans et al. 2013) because of the Earth occultations. Here, we are 
able to establish this function using the whole HST orbit and con¬ 
firm that it is well approximated by a 4th-order polynomial. The 
logarithmic ramp is necessary to reproduce the increase in flux seen 
by visual inspection of the time-series. 

This basic detrending reveals correlated structures in the pho¬ 
tometry suggesting an imperfect correction. It is indeed unlikely 
that 26-hr of HST continuous observations can be accurately mod¬ 
elled with a single baseline model (constrained by fixed coeffi¬ 
cients) as simple as the one described above (see Figure 5). Since 
we do not know a priori if a transit exists in the data and what its 
location could be, we analyse the time-series by separating them 
in different segments spanning four to eight orbits, similar to other 
HST exoplanet programs. All segments are analysed in the same 
MCMC fit, with their own baseline model coefficients, but retain¬ 
ing the same functional form (orbital phase dependent polynomial 
and linear trend). We first split our photometric time-series into two 
segments. The logarithmic ramp term is needed for the first seg¬ 
ment only. We then split our data into three segments. Since the 
ramp seen on the whole photometry affects the first 0.6 days of 
our data, the baseline models for the first two segments include a 
logarithmic ramp model. The time-series are finally split into four 
segments. For this 4-segment MCMC fit, only the first two files 
have a baseline including a ramp model. The results for these 3 
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MCMC analyses are shown in Figures 6, 7 and 8. We find that 
splitting the photometry in segments significantly decreases the 
Bayesian Information Criterion (BIC, Schwarz 1978) as compared 
to a unique fit to the whole photometry, despite the significant in¬ 
crease in free parameters. The behaviour of HST/STIS combined 
systematics slightly evolve from one orbit to another, making it 
difficult to model 16 orbits with a unique set of baseline coeffi¬ 
cients. As expected, most structures found in the full photometry 
fit have vanished, but one, centred on To = 2456482.195 BJD 
with a ~0.15-day duration (first-to-fourth contact). This transit-like 
structure has a location, depth and duration that remains consistent 
whether the data are unsegmented or split in 2, 3 or 4 segments (see 
Figures 5,6, 7 and 8). We notice another structure in the time-series 
with a 0.2 to 3-er significance (depending on the number of seg¬ 
ments) and located at To = 2456481.9 BJD. However, we find the 
duration and depth of this structure to change significantly with the 
baseline model and the number of segments used in the analysis, 
which suggests an instrumental origin. For comparison purposes, 
the transit-like pattern centred on To = 2456482.195 BJD is de¬ 
tected at the 7.5 -a level. We show in Table 1 the MCMC fit results 
obtained for each analysis (unsegmented and segmented by 2, 3 and 
4 parts). 

We then investigate whether the inclusion of ancillary param¬ 
eters measured during the data reduction helps in improving the fit. 
We find that including the median overscan value as a linear coeffi¬ 
cient in the baseline function improves the BIC by 3%. The inclu¬ 
sion of linear combinations of the spectrum drifts in the dispersion 
and cross-dispersion directions increases the BIC by 5%. Our final 
baseline model therefore consists of HST’s orbital phase’s fourth 
order polynomial, a time- and overscan-dependent linear trends. 
Adding higher-order terms to these functions does not improve the 
BIC. 

We find a consistent exposure-to-exposure photometric preci¬ 
sion of 110 to 120 ppm depending on the segment, significantly 
above the Poisson-limited precision of 20 ppm per exposure. The 
origin of this large discrepancy is not known but could be due to 
a larger than expected high-frequency instrumental noise due to 
shutter-timing jitter or could also have a contribution from stellar 
granulation or pulsations. We measure low contamination from cor¬ 
related noise (Pont et al. 2006) that we estimate between 22% and 
34% following Gillon et al. (2010) depending on the segment. Fig¬ 
ure 9 shows the actual decrease of the photometric RMS with the 
time-series bin size (black points). The theoretical photon noise is 
shown as a solid red line. We also indicate the expected decrease 
in Poisson noise normalised to the individual exposure time (in red 
dotted line). Despite the larger discrepancy between Poisson noise 
and our actual precision, Figure 9 shows that this extra noise com¬ 
ponent has mostly a white behaviour that does not affect our sci¬ 
ence goal (transit detection). The improvement in precision is even 
accentuated for the larger bin sizes (60 min and more), meaning 
that a large part of the extra noise noticed at higher frequencies 
reduces significantly for larger time bins, which correspond to typ¬ 
ical transit durations. We adjust our individual error bars taking 
into account these extra white and red noise components within the 
MCMC framework. 

In the following we focus on the observed transit structure 
centred on 2456482.2 BJD. The second half of the transit structure 
benefits from HST continuous monitoring with no interruptions due 
to the SAA, which enabled the continuous monitoring of the tran¬ 
sit egress. The first half of the transit window suffers from regular 
interruptions, which prevented us from monitoring the ingress. We 
perform a series of fits with different functional forms of the base¬ 


line models, including orbital phase polynomial orders from 4 to 7, 
overscan polynomial orders from 0 to 5 as well as the removal of 
the linear trend. The resulting transit shape (depth and duration) is 
found to be consistent at the l-cr level in each case. The adopted 
baseline model (in the BIC sense) remains the same as the one 
used for the global fit described above. We find no correlation of 
the transit signature in the photometry with any ancillary param¬ 
eter (spectrum position, background or overscan values measured 
during data reduction). A zoom of the transit pattern is shown in 
Figure 10. 

2.2 July 2014 dataset 

The transit pattern detected in our July 2013 dataset prompted us 
to request new observations to confirm the transit signal, based on 
the radial-velocity period of 3.2357 days (Dumusque et al. 2012) 
and the transit centre computed from our STIS July 2013 observa¬ 
tions. Nine orbits of HST DD time (PID 13927) were awarded to 
re-observe a Cen B for 13.5 hours from 11:30 UT on 28 July 2014 
to 01:35 UT on 29 July 2014. 

The data reduction and photometric analysis employed for this 
dataset are identical to the July 2013 observations. The raw pho¬ 
tometry is shown in Figure 11. We find that the first orbit of this 
visit exhibits a pronounced trend, likely due to the telescope set¬ 
tling on a new attitude. We do not notice such behaviour for the July 
2013 dataset. We thus elect to discard the first 56 min (first orbit) 
of this visit. Ignoring the first orbit of a visit is a common proce¬ 
dure but is not found to be necessary for all STIS observations. The 
calibrated time-series for this second visit are shown in Figure 12. 
We measure an RMS of 116 ppm and a red noise contribution of 
20%, in excellent agreement with the first visit’s photometric prop¬ 
erties, emphasising on the consistency of saturated spectroscopy 
with STIS. 

We do not find any structure matching the transit pattern de¬ 
tected in the first dataset, while the photometric precision achieved 
for these 9 orbits would have enabled a clear detection of a 3.8-hr 
long, 100-ppm deep transit-like signature in the data at more than 
the 5 -<t level. 


3 NO TRANSIT FOR ALPHA CEN B B 

In total, we observed a Cen B for almost 40 hours in a consistent 
way with HST/STIS. Both the July 2013 and 2014 datasets exhibit 
low-levels of correlated noise and a 6-s RMS of ~ 115 ppm allow¬ 
ing us to confidently determine the a posteriori transit probability 
of the Earth-mass exoplanet a Cen B b. We combine in the fol¬ 
lowing both the 2013 and 2014 photometric datasets in an MCMC 
framework with the aim of determining to which extent our dataset 
rules out a transit of a Cen B b with the published parameters of 
Dumusque et al. (2012). 

We include both the published radial-velocity data and the 
STIS photometry described above as input data to the MCMC 
fit. We model the radial velocities following the appendix of 
Dumusque et al. (2012). We use the same functional form for the 
detrending model including the rotational activity, magnetic cy¬ 
cles and binary motion. The detrended radial-velocity datapoints 
are included in the MCMC fit and their errors include the model 
coefficient uncertainties. For the transit model, we compute the 
quadratic limb-darkening coefficients from a Kurucz stellar atmo¬ 
sphere model (Kurucz 1979). We first adopt a 5,250K, log</=4.5 
and [Fe/H]=0.2 model, which represents the closest match to 
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Figure 5. Detrended photometry - unsegmented. Data are binned per 45 minutes. The best fit model is superimposed in red. 



Figure 6. Detrended photometry - 2 segments. Data are binned per 45 minutes. The best fit model is superimposed in red. 
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Figure 7. Detrended photometry - 3 segments. Data are binned per 45 minutes. The best fit model is superimposed in red. 


! T= =T_T_ T! 

t T X_ T _T-r_T_T 


. 1 . 1 . 1 .. 

b 

b 

W 

b 

b 

' H 

b 

W 

b 

b 

b 

M 

r wj 


1.0003 
1.0002 
1.0001 
1.0000 
0.9999 
0.9998 
0.9997 
0.9996 

6481.3 6481.4 6481.5 6481.6 6481.7 6481.8 6481.9 6482.0 6482.1 6482.2 6482.3 

Time [BJD-2450000] 




- T T I T I ... 


±± — I — X - * 

T X 1 -*■ 

£ * I I 1 1 



Figure 8. Detrended photometry - 4 segments. Data are binned per 45 minutes. The best fit model is superimposed in red. 
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Table 1. Transit candidate fit results. MCMC fit results obtained for each analysis (unsegmented and segmented by 2, 3 and 4 parts). Values are the median 
of the posterior distributions for each perturbed parameter indicated on the left column, along with the 1 -<j credible intervals. BIC values are calculated for the 
combined baseline+transit model. 



Figure 9. RMS of binned residuals versus bin-size. The first panel corresponds to the first segment, the fourth panel to the last segment. The theoretical 
photon noise is shown as a solid red line. The photometric RMS is indicated vs. the time-series bin size (black points). The expected decrease in Poisson noise 
normalised to an individual exposure (6s) is shown as a red dotted line. 


a Cen B's properties. We use the model output file to calculate 
the weighted averaged centre-to-limb intensity profile using STIS 
G750M’s 6094 setting wavelength coverage. We then fit the profile 
with a two-parameter quadratic limb-darkening law (Claret 2000). 
We repeat the same procedure for the 6 neighbouring stellar atmo¬ 
sphere models in the 3-dimensional grid and use the uncertainties 
on a Cen B's stellar properties to estimate individual errors on the 
limb-darkening parameters. We find the quadratic coefficients to be 
ui = 0.504 ± 0.010 and u 2 = 0.236 ± 0.008. The MCMC fit has 
the following jump parameters: the time of RV maximum To, the 
orbital period P, the parameter K' = Ky/1 — e 2 P 1 ^ 3 , where K 
is the radial-velocity semi-amplitude, the two parameters fe cos w 
and x/esinu;, the planet/star radius ratio and the impact 

parameter b. We also add the linear combinations ci = 2«i + u 2 
and c 2 = ui — 2u 2 (where u\ and u 2 are the quadratic coeffi¬ 
cients determined above) as jump parameters to minimise the cor¬ 
relations of the resulting uncertainties (Holman et al. 2006). We im¬ 
pose Gaussian priors on the limb-darkening coefficients. 

We note that neglecting the 2014 dataset yields an orbital 
solution providing a satisfying fit to both the spectroscopic and 
photometric data, provided the orbit has a significant eccentricity 
(e = 0.54) and a longitude of periastron u> ~ 270 degrees. This 
orbital configuration allows e sin to to be significantly offset from 
0, thus providing leverage to the fit in increasing the transit du¬ 
ration while leaving the transit To almost unchanged. If this tran¬ 
sit was caused by a Cen B b, the transit depth would translate to 
a planetary size of 0.92 ± 0.06 Earth radii, which, combined to 
our revised mass measurement of 1.01 ± 0.09 Earth masses, re¬ 
sults in a density that is 30% larger than the Earth, similar to that 
of Mercury (Haucketal. 2013). Using specific dissipation func¬ 
tion values of Goldreich & Soter (1966) to estimate circularisation 
timescales, we estimate eccentricity damping times of less than 
100 Myr for a planet on such a large eccentricity orbit, which is 
significantly shorter than the age of 5-6 Gyr of the a Cen system 
(Eggenberger et al. 2004). While not impossible, it is unlikely that 


the eccentricity of such a small planet would be maintained at ec¬ 
centricities > 0.5 as suggested by the transit signal duration. Com¬ 
bining both HST datasets rules out the link between the transit-like 
pattern seen in July 2013 and a Cen B b published orbital param¬ 
eters. Note however that the 2014 dataset probes the repeatability 
of the transit pattern within a 2-a window based on the published 
orbital period value. 

Our MCMC analysis also allows us to compute the a pos¬ 
teriori transit probability of a Cen B b based on the parameters 
of Dumusque et al. (2012). In other words we estimate to which 
level our combined HST dataset discards the transit of a Cen B b 
as published. Using the published orbital parameters, the geomet¬ 
ric a priori transit probability of a Cen B b assuming a circular 
orbit is 9.5%. For the purpose of this analysis, we put a normal 
prior on the planetary radius N(0.97,0.10 2 ), which is motivated 
by the ~1 Earth-mass of the planet. We also include the radial- 
velocity data and consider only circular orbits. We find that for 
3.4% of the MCMC accepted steps, a transiting configuration is 
found requiring a To ~ 2455281.55 and P ~ 3.2375, at the 
edge of our observing window. These values have to be com¬ 
pared to the transit ephemerides of To = 2455281.24 ± 0.17 and 
P = 3.2357 ± 0.0008 computed from the radial-velocities alone, 
prior to the HST observations. The a posteriori transit probability 
is thus 0.034 x 0.095 = 0.32%. Our observations thus allows to 
reduce the transit probability of a Cen B b by a factor 30. 


4 POSSIBLE ORIGIN OF THE 2013 TRANSIT PATTERN 

We explore in the following the possibility that the July 2013 tran¬ 
sit pattern is due to stellar variability, instrumental systematics or 
cause by a background eclipsing binary. We do not find any temper¬ 
ature or HST orbital dependent parameter, nor X/Y spectral drifts 
to correlate with the transit pattern. The transit candidate duration 
of 3.8 hours is 2.4 times longer than the HST orbital period, mak- 
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Figure 10. Corrected HST/STIS corrected light-curve centered on the transit signal. Data are binned per 45 minutes. Unbinned datapoints are shown in 
grey. Our best-fit MCMC model is superimposed in red. 


ing the transit pattern unlikely to be attributable to HST instrumen¬ 
tal systematics. As the detector is consistently saturated during all 
of our observations, we also find it unlikely that saturation is the 
origin of the transit signal. Another possible explanation could be 
stellar variability. However, the duration of the transit candidate 
(3.8-hr) is not consistent with the stellar rotational period of 36.2 
days (DeWarf et al. 2010), to enable a spot (or group of spots) to 
come in and out of view. In such a case, star spots would change the 
overall observed flux level and produce transit-shape signals, as is 
the case for stars having fast rotational periods (e.g. Harrison et al. 
2012). We finally find no sign of light contamination from the A 


component, which we specifically avoided during the preparation 
of both 2013 and 2014 visits. Our data reduction and analysis is 
straightforward, including only basic detrending variables related 
to the HST orbital phase and spectrum position. We find the transit 
pattern to be robust to the detrending approaches we employ. 

We then estimate the probability that the observed transit of 
91 ppm is due to a background eclipsing binary (BEB). We first 
estimate up to which magnitude a diluted BEB can cause this signal 
by applying 


Sbeb =Seb■ 


Feb 


FeB + Ffarget 


Am 


= Seb • 2.512' 


( 1 ) 
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Figure 11. July 2014 Raw HST/STIS photometry. The single-baseline best-fit MCMC model is superimposed in red. 



Time [BJD-2450000] 

Figure 12. July 2014 Detrended photometry - 2 segments. Data are binned per 45 minutes. See text for details. 


whereby we express the flux ratio via the magnitude difference Am 
between the BEB and the target star. We now allow all 8eb £= (0,1] 
and find that the target star with V = 1.33 can dilute BEBs of up 
to V < 11.5. Fainter BEBs cannot cause a deep enough signal, and 
only faint BEBs with significant transit depths or bright BEBs with 
grazing eclipses can cause the required depth. Secondary eclipses 
will further drastically reduce the fraction of reasonable systems. 
However, here we allow every scenario with V < 11.5 and thus 
significantly overestimate the total number of possible BEBs. Us¬ 


ing the TRILEGAL (Girardi et al. 2005) galaxy model we find 
a distribution of around 1000 star systems with V < 11.5 in 
1 sq. deg. around the field of view, leading to 0.8% probability to 
randomly have any of these in the 52 x 2 sq. arcsec slit. The prob¬ 
ability to have any kind of detached EB among all star systems is 
roughly 1%, shown by our simulations and confirmed by the results 
from Kepler (Slawson et al. 2011). Thereby, we again overestimate 
the risk of such BEBs, since triple and higher-order systems ad¬ 
ditionally dilute any eclipse such that the magnitude limit for the 
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BEB would be lower, leaving significantly less systems. The prob¬ 
ability for the signal to be caused by a BEB is hence lower than 
0.008%. 

A possible remaining explanation for the observed signal is 
thus a planetary transit whose orbital period is significantly dif¬ 
ferent from a Cen B b’s, as published in the discovery paper. We 
use our HST transit detection to constrain the orbital period of 
this putative planet and the observing window where another tran¬ 
sit would happen, which would confirm the repeatability of the 
signal. The host star, Alpha Cen B, benefits from extensive char¬ 
acterisation. Both the stellar mass (M* = 0.934 ± O.OO6M 0 ) 
and radius (7?* = 0.862 ± 0.0057 ?q) are known to a remark¬ 
able precision (Kervella et al. 2003). We thus know the stellar 
density, which can also be constrained by the transit shape (e.g. 
Seager & Mallen-Ornelas 2003, assuming a circular orbit.). Com¬ 
bining our knowledge from the host star and the 2013 transit 
shape allows us to constrain the orbital period, impact parame¬ 
ter and eccentricity, albeit in a degenerate way. We employ our 
MCMC framework described above to explore what orbital peri¬ 
ods/eccentricities/impact parameters are allowed by the HST 2013 
transit lightcurve alone. We find that the STIS photometry yields an 
orbital period 2-a upper limit of 20.4 days (the median of the poste¬ 
rior being 12.4 days), with a modest impact parameter ~ 0.0 — 0.3 
and an eccentricity 2 -a upper limit of 0.24. 

The transit pattern detected in our HST/STIS data motivates 
further photometric follow-up observations of a Cen B to con¬ 
firm its repeatability. It is unlikely that the mass of an Earth-size 
planet orbiting a Cen B with a significantly longer orbital period 
than a Cen B b could be constrained in the near future. Since the 
radial-velocity semi-amplitude decreases with 1/P 1 / 3 , an Earth- 
mass planet on a e.g. 15-day orbit would yield a radial-velocity 
amplitude of 0.3 m/s, compared to K=0.5 m/s in the case of the 
3.24-day a Cen Bb. The photometric follow-up is equally compli¬ 
cated. Ground-based facilities do not have the ability to 1) con¬ 
tinuously monitor the 2-3 weeks window and 2) to reach the 30 
ppm per 3.8 hour precision necessary to detect a 90-ppm transit 
(assuming an Earth-like interior) at the 3-cr level. The MOST satel¬ 
lite (Walker et al. 2003) performed several observations of bright 
stars. Considering the a Cen system, the transit signature would be 
diluted by a factor of three, in the combined light of the close A 
and B components, thus rendering its confirmation extremely diffi¬ 
cult. a Cen B’s flux in Spitzer ’s IRAC channels 1 and 2 is ~371 and 
~220 Jy respectively, more than ten times IRAC’s saturation levels. 
Our observing strategy and data reduction technique demonstrates 
that HST/STIS is able to detect 100-ppm transits at 5-<r or more 
for stars that exceed the saturation level of the detector by a factor 
5. The characterisation of small exoplanets orbiting stars as bright 
as a Cen B will continue to pose significant challenges in the fu¬ 
ture unless next-generation of telescopes have the ability to employ 
successful observing strategies for highly-saturated stars. 


ACKNOWLEDGMENTS 

We are indebted to Shelly Meyett and Tricia Royle for their assis¬ 
tance in the planning and executing of these observations. We thank 
Tom Ayres for his help in preparing our phase 2 program, as well as 
Adrian Barker, Samantha Thompson and Julien de Wit for discus¬ 
sions. We thank the anonymous referee for a detailed and helpful 
review that improved the paper. We are grateful to STScI Direc¬ 
tor Matt Mountain for awarding Director’s Discretionary Time for 
both programs. Based on observations made with the NASA/ESA 


Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program PID13180 and 
13927. Support for this program was provided by NASA through 
a grant from the Space Telescope Science Institute, which is oper¬ 
ated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. This work has been car¬ 
ried out within the frame of the National Center for Competence in 
Research PlanetS supported by the Swiss National Science Funda- 
tion (SNSF). D.E., C.L., F.P., and S.U. acknowledge the financial 
support of the SNSF. A. T. is a Swiss National Science Foundation 
fellow under grant number P300P2-147773. 


REFERENCES 

Batalha N. M., Rowe J. E, Bryson S. T., Barclay T., Burke C. J., Caldwell 
D. A., Christiansen J. L., 2013, ApJS, 204, 24 
Boisse I.. Bouchy F., Hebrard G., Bonfils X., Santos N., Vauclair S., 2010, 
European Planetary Science Congress 2010, p. 151 
Claret A., 2000, A&A, 363, 1081 

Demory B.-O., Gillon M.. Deming D., Valencia D., Seager S., Benneke B., 
Lovis C., Cubillos P., Harrington J., Stevenson K. B., Mayor M., Pepe F., 
Queloz D., Segransan D., Udry S., 2011, A&A, 533, 114 
DeWarf L. E„ Datin K. M„ Guinan E. F„ 2010, ApJ, 722, 343 
Dragomir D., Matthews J. M., Eastman J. D., Cameron C., Howard A. W., 
Guenther D. B., Kuschnig R., Moffat A. F. J., Rowe J. F., Rucinski S. M., 
Sasselov D., Weiss W. W., 2013, ArXiv e-prints 
Dumusque X., Pepe F., Lovis C., Segransan D., Sahlmann J., Benz W., 
Bouchy F., Mayor M., Queloz D., Santos N., Udry S., 2012, Nat, 491, 
207 

Eastman J., Siverd R., Gaudi B. S., 2010, PASP, 122, 935 
Eggenberger P., Charbonnel C., Talon S., Meynet G., Maeder A., Carrier 
F„ Bourban G„ 2004, A&A, 417, 235 
Evans T. M., Pont F., Sing D. K.. Aigrain S., Barstow J. K., Desert J.-M., 
Gibson N., Heng K., Knutson H. A., Lecavelier des Etangs A., 2013, 
ApJ, 772, L16 

Gilliland R. L., 1999, Technical report. Instrument Science Report STIS 
1999-05 (Baltimore, MD: STScI) 

Gilliland R. L„ Goudfrooij P., Kimble R. A., 1999, PASP, 111, 1009 
Gillon M., Bonfils X., Demory B.-O., Seager S., Deming D., Triaud 
A. H. M. J., 2011, A&A, 525, A32 

Gillon M., Lanotte A. A., Barman T., Miller N„ Demory B.-O., Deleuil 
M., Montalban J., Bouchy F., Collier Cameron A., Deeg H. J., Fortney 
J. J., Fridlund M., Harrington J., Magain R, Moutou C., Queloz D., Rauer 
H., Rouan D., Schneider J., 2010, A&A, 511,3 
Gillon M., Pont F., Demory B.-O., Mallmann F., Mayor M., Mazeh T., 
Queloz D., Shporer A., Udry S., Vuissoz C., 2007, A&A, 472, L13 
Gillon M., Triaud A. H. M. J., Fortney J. J., Demory B.-O., Jehin E., Lendl 
M., Magain R, Kabath R, Queloz D., Alonso R., Anderson D. R., Collier 
Cameron A., Fumel A., Hebb L., Hellier C., Lanotte A., Maxted P. F. L., 
Mowlavi N„ Smalley B„ 2012, A&A, 542, A4 
Girardi L., Groenewegen M. A. T., Hatziminaoglou E., da Costa L., 2005, 
A&A, 436, 895 

Goldreich P., Soter S., 1966, Icarus, 5, 375 

Harrison T. E., Coughlin J. L., Ule N. M., Lopez-Morales M., 2012, AJ, 
143,4 

Hatzes A. R, 2013, ApJ, 770, 133 

Hauck S. A., Margot J.-L., Solomon S. C., Phillips R. J., Johnson C. L., 
Lemoine F. G., Mazarico E., McCoy T. J., Padovan S., Peale S. J., Perry 
M. E., Smith D. E., Zuber M. T., 2013, Journal of Geophysical Research 
(Planets), 118, 1204 

Heng K., Showman A. R, 2014, ArXiv e-prints 

Holman M. J., Winn J. N., Latham D. W., O’Donovan F. T., Charbonneau 
D„ Bakos G. A., Esquerdo G. A., Hergenrother C., Everett M. E., Pal A., 
2006, ApJ, 652, 1715 


HST search for the transit of a Cen Bb 11 


Howard A. W., Marcy G. W., Johnson J. A., Fischer D. A., Wright J. T., 
Isaacson H., Valenti J. A., Anderson J., Lin D. N. C., Ida S., 2010, Sci¬ 
ence, 330, 653 

Huitson C. M., Sing D. K., Vidal-Madjar A., Ballester G. E., Lecavelier 
des Etangs A., Desert J.-M., Pont F., 2012, MNRAS, p. 2826 
Kervella R, Th Venin F., S Gransan D., Berthomieu G., Lopez B., Morel 
R, Provost J., 2003, A&A, 404, 1087 
Knutson H. A., Charbonneau D., Allen L. E., Burrows A., Megeath S. T., 
2008, ApJ, 673, 526 
Kurucz R. L., 1979, ApJS, 40, 1 

Mayor M., Marmier M., Lovis C., Udry S., Segransan D., Pepe F., Benz 
W., Bertaux J. ., Bouchy F., Dumusque X., Lo Curto G., Mordasini C., 
Queloz D., Santos N. C., 2011, ArXiv e-prints 
Pont F., Zucker S., Queloz D., 2006, Monthly Notices of the Royal Astro¬ 
nomical Society, 373, 231 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P, 1992, Nu¬ 
merical recipes in FORTRAN. The art of scientific computing 
Schwarz G., 1978, The Annals of Statistics, 6, 461 
Seager S., Deming D., 2010, ARA&A, 48, 631 
Seager S., Mallen-Omelas G., 2003, ApJ, 585, 1038 
Sing D. K., Desert J.-M., Lecavelier Des Etangs A., Ballester G. E., Vidal- 
Madjar A., Parmentier V., Hebrard G., Henry G. W., 2009, A&A, 505, 
891 

Slawson R. W., Prsa A., Welsh W. F., Orosz J. A., Rucker M., Batalha N., 
Doyle L. R., Engle S. G., 2011, AJ, 142, 160 
Walker G., Matthews J., Kuschnig R., Johnson R., Rucinski S., Pazder J., 
Burley G., Walker A., Skaret K., Zee R., Grocott S., CaiToll K., Sinclair 
P, Sturgeon D., Harron J., 2003, PASP, 115, 1023 
Winn J. N., Matthews J. M., Dawson R. I., Fabrycky D., Holman M. J., 
Kallinger T., Kuschnig R., Sasselov D., Dragomir D., Guenther D. B., 
Moffat A. F. J., Rowe J. F., Rucinski S., Weiss W. W., 2011, ApJ, 737, 
L18 


This paper has been typeset from a TgX/ ET^X file prepared by the 
author. 


